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We show that chiral symmetry and gauge invariance enforce relations between the short-distance 
physics that occurs in a number of electroweak and pionic reactions on hght nuclei. Within chiral 
perturbation theory this is manifested via the appearance of the same axial isovector two-body 
contact term in ■K~d — > nwy, p-wave pion production in A^A'' collisions, tritium f3 decay, pp fusion, i/d 
scattering, and the hep reaction. Using a Gamow- Teller matrix element obtained from calculations 
of pp fusion as input we compute the neutron spectrum obtained in vr^d — > nwy. With the short- 
distance physics in this process controlled from pp — > de'^Ue the theoretical uncertainty in the nn 
scattering length extracted from n~d — > nn-y is reduced by a factor larger than three, to < 0.05 fm. 
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In Quantum Chromodynamics (QCD) the up and 
down quark masses are much smaller than the typical 
scale of hadron masses, Mqcd ^ 1 GeV/c^, and so the 
QCD Lagrangian has an approximate S'[/(2)l x 5J7(2)r 
symmetry. This symmetry is spontaneously broken by 
the QCD vacuum, leaving only 5'[/(2)isospin as a man- 
ifest symmetry of low-energy hadronic processes. The 
spontaneously-broken symmetry results in the existence 
of three (pseudo-)Nambu-Goldstone bosons, the pions, 
and means that the quantum-field-theoretic operator cor- 
responding to the Noether currents of the broken gen- 
erators changes pion number by one. External elec- 
troweak fields couple to these currents, and so the chiral 
SU{2)i^ X SU{2)b. symmetry of QCD imposes relation- 
ships between the pion-baryon couplings and the axial 
couplings of baryons measured in electroweak processes. 
One of these is the Goldberger-Treiman relation: 
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which expresses the nNN coupling constant g-^NN in 
terms of the nucleon's axial coupling gp^ and the pion 
decay constant f-^. [M is the nucleon mass.) 

The symmetries of QCD have other consequences for 
low-energy reactions involving pions. Introducing a pho- 
ton field into the theory according to the dictates of 
f7(l)om gauge invariance shows that the leading piece of 
the amplitude for the photoproduction of 7r+ or 7r~ at 
threshold is given by the KroU-Ruderman (KR) term: 
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In the singie-nucleon sector the combination of J7(l)cm 
and QCD's spontaneously-broken SU (2)l x SU (2)r en- 
forces relations like Eqs. and In this paper we 
explore the analog of these relations in the two-nucleon 
sector. We first exhibit a relationship between the short- 
distance physics in NN NNtt reactions and elec- 



troweak processes such as pp fusion, tritium beta de- 
cay, and muon absorption on deuterium fi~d — > nnv^. 
We then show that the same two-nucleon physics oc- 
curs in some reactions involving photons, in particular 
in the process ■K~d —f nn^. The two-nucleon coupling 
in this reaction is then directly proportional to the two- 
nucleon coupling of the axial current. Knowledge of one 
coupling therefore yields the other, so that precise in- 
formation about n~d — ^ nnj (or its crossed partner 
jd — > nnTT"*") constrains pp fusion. Conversely, calcu- 
lations of tritium beta decay and pp fusion contribute to 
model-independent predictions for n^d nn"/. 

Equations (Q) and ^ are exact only in the chiral limit 
of QCD, i.e., when to„ = = 0. In our world they are 
satisfied to within 2% and 6% respectively, since they 
have corrections proportional to m„ and m^. The tool 
for systematically computing such corrections is chiral 
perturbation theory (xPT) (see Ref. Q] for a review). 
XPT is an effective field theory (EFT) in which ampli- 



tudes are obtained in an expansion in x 
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Since 



its initial development 0, Q xPT has been successfully 
applied to a number of processes in the meson-meson and 
meson-nucleon systems. But few-nucleon systems involve 
states which are bound by the strong interaction, and so 
the xPT expansion must be modified. Weinberg 0| has 
pointed out that the xPT expansion can be applied to the 
operator O used in the computation of a matrix element 
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for reactions involving such systems. The result — 
rovided consistently-computed wave functions {ipij) [M 
are used — is an expansion of M in powers of x that 
incorporates the model-independent consequences of the 
(broken) chiral symmetry of QCD into calculations of 
processes in few-nucleon systems. 

The operator O in Q is computed up to a given order 
in the xPT expansion in powers of x, using the Feynman 
rules derived from the xPT Lagrangian. At leading-order 
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the relevant terms for one nuclcon are: 

Cn =N^[iv-D + gAS-u]N (4) 

where v is the nucleon four- velocity, D is a, (chiral) covari- 
ant derivative. S is the Pauli-Lubanski spin vector, and 
5a parameterizes the unknown short-distance physics of 
the nucleon's axial-current matrix element. Here is an 
axial four-vector which contains the pion field, and when 
expanded takes the form: 

where (^m) external vector (axial) field. The 

Goldberger-Treiman and KR relations 1^ and (O follow 
directly from inserting jSJl into the Lagrangian Q and 
examining the relevant ttNN and jttNN couplings. 

The axial four-vector also appears when pions, pho- 
tons, and axial fields couple to pairs of nuclcons. The 
leading effects of this can be represented by one term in 
the chiral Lagrangian 0, |^ : 

Cnn = -2diN''S ■uNN^'N. (6) 

In principle the low-energy constant (LEG) di can be 
evaluated by lattice methods [3- As with most single- 
nucleon sector LECs, however, its value has at present 
only been obtained from experimental data. Once the 
value of the LEG di has been extracted from one pro- 
cess we can use the result to predict other observables 
on which it has an impact. Equation © thus encodes 
model-independent correlations between diiferent reac- 
tions which share the same short-distance physics be- 
cause of the symmetries of the chiral Lagrangian. For the 
reactions discussed here the short-distance NN physics 
is associated with '^S'l ^Sq transitions. It is convenient 
to parameterize it by a dimensionless constant, d, which 
is obtained from di, and from which the scale Mf^ has 
been removed 0, E| ■ In two-nucleon systems d plays a 
role analogous to that of gA in the single- nucleon sector. 

Model-independent correlations involving d have been 
used to predict the rates of reactions which are important 
for the production and detection of solar neutrinos. In 
Ref. the constant d was fixed using the well-measured 
tritium beta-decay half-life. This allowed the authors 
of that paper to predict the rate of two reactions which 
are important members of the chain that produce en- 
ergy in main-sequence stars and result in the emission of 
solar neutrinos: the proton-fusion process pp de^v^. 
and the hep process '^Hep — > '^Hce'^i'e- Importantly, d 
also parameterizes the short-distance physics in the v(p)d 
breakup reactions ^3 that facilitate the flavor decompo- 
sition of the solar-neutrino signal from data obtained at 
the Sudbury Neutrino Observatory (SNO) [13 ■ Since the 
pp fusion, hep, and v{v)d breakup reactions occur at too 
low a rate to be reliably measured in the laboratory, the 
correlations exhibited in Refs. 0.llll| provide important, 



model-independent information on them. (This connec- 
tion has also been explored in potential model calcula- 
tions [ill and pionless EFT The same Gamow- 
Teller '^Si ~* ^Sq transition occurs in /i^d nnv^, and 
the LEG d plays a role there too 0| . 

If we now insert the expression jSJ into the Lagrangian 
© we see that if a pion is produced via the reaction 
NN NNtt, and is in a p-wave relative to the NN 
system, then the LEG d parameterizes the short-distance 
physics of the process. This is the two-body analog of 
the Goldberger-Treiman relation and it establishes 
an intriguing connection between the reactions involved 
in the production and detection of solar neutrinos and 
the physics of pion production in NN collisions. 




FIG. 1: The piece of the leading chiral 3NF which involves d. 

In systems with A > 2, d plays a subtle, but important, 
role in nuclear binding through its generation of a three- 
nucleon force (3NF). As pointed out in Ref. 0|, a 3NF is 
obtained if a virtual pion is produced via © and then ab- 
sorbed on a third nucleon via the leading Lagrangian Q). 
The resulting diagram is a leading non- vanishing term in 
the chiral expansion of the 3NF [3 ^| , and is depicted 
in Fig. HI D etermining d accurately from 3N data may be 
difficult [l7| J but the correlations mandated by the chiral 
symmetry of QGD mean that we can use the values of d 
extracted from electroweak processes — e.g., in Ref. — 
as input to fix the size of the 3NF depicted in Fig.^ This 
has the advantage that the kinematics at which the vir- 
tual pion is emitted in Fig. ^ (pion energy w 0) matches 
the kinematics at which reactions like pp — > de^v^ oc- 
cur quite well — certainly much better than it does the 
kinematics of NN NNtt, which was used to constrain 
the 3NF in Ref. [3- This should reduce the impact that 
higher-order xPT corrections have on predictions for this 
piece of the 3NF. Qualitative comparison of the values 
extracted for d in Refs. and ^3 is encouraging, but 
for a quantitative comparison tritium /3-decay and three- 
nucleon bound state and scattering calculations will have 
to be done with consistent wave functions and regulators. 

Here, we are less ambitious. We will instead exploit 
another correlation that results from substitution of Q 
into — one involving photons. We examine the reac- 
tion TT^d nwy, which has played a significant role in 
determinations of the nn scattering length (a„„). Data 
from radiative pion capture experiments dominates the 

accepted value a„„ ~ — 18.59±0.40 fm [II El SHI [13. 

Final-state interactions in nd breakup experiments can 
also be used to measure a„„. While one such recent 
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determination is in agreement with the accepted value, 
another measurement of a„„ usin g th e same reaction dis- 
agrees with it by more than 4cr [23 • We will not com- 
ment on this further here, but will focus on showing that 
there is a correlation between 7r~c? nnj and pp fusion. 
This correlation, together with knowledge of the pp fu- 
sion rate, reduces the impact of short-distance physics on 
the value of a„„ extracted from Tr~d —^ nn'y data. This is 
important because previous extractions of a„„ from this 
process have shown that ±0.3 fm out of the to- 

tal uncertainty of ±0.4 fm in the final result comes from 
the short-distance piece of the NN wave function. Up 
until now it was not clear how to remove this source of 
uncertainty, and so reduce the a„„ error below ±0.3 fm. 

We use wave functions that are derived in a man- 
ner consistent with chiral symmetry. For the deuteron 
S'-state, we start from the asymptotic wave function 
M(o)(r) = Ase-''^ where As = 0.8846 fm^^/^ ^j^e 
asymptotic normalization and 7 = \/MB = 45.7 MeV 
the binding momentum |24l |. For low-energy nn scatter- 
ing the starting point is asymptotic wave functions with 
phase shifts derived from the effective-range expansion, 
i.e., v^^^{r) = sin(pr + S)/smS, where p is the relative 
nn cm. momentum and pcotS = — l/ftnn + j'^oP^ (''0 is 
the effective range) . In both cases the wave functions are 
calculated from r = 00 down to a matching point R us- 
ing the one-pion-exchange potential. For r < R we then 
assume a spherical well potential and match to the r > R 
wave function at r = i?. Further details can be found in 
Refs. jiEEEIi^- The pp state needed for the calculation 
of pp — > de'^i'e is calculated similarly, with the asymp- 
totic wave function given by zero-energy Coulomb wave 
functions, and the pp scattering length with respect to 
the Coulomb potential, ag, = -7.8196 ± 0.0026 fm 
as input. Higher-order electromagnetic corrections can 
be included, and give a < 1% effect 0|, but for our 
present purposes they are not important, since they alter 
only the long-distance properties of the pp wave function. 

The matrix element for proton fusion is given by |2^ : 
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up to corrections of relative order x^- In Eq- u{r) and 
i;^(r) are the deuteron S'-state and pp ^Sq zero-energy 
scattering (including Coulomb) wave functions. The ma- 
trix element for the KR contribution to Tr~d nnj in 
the final-state interaction (FSI) region (see Fig. ^ is 
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FSI 



C / dru(r)jo 



v{r), 



(8) 



where C is a (known) constant, k the photon cm. mo- 
mentum, jo a spherical Bessel function, and v{r) the full 
(energy-depend ent) nn ^Sq scattering wave function in 
the notation of |22| . Equation JHJ is also of accuracy ■ 
Figure [3 then shows that the Gamow- Teller matrix 
element O is correlated with the FSI peak height in 



TT~d nn7, which here, and in Fig. O below, we have 
computed using the 0{x'^) xPT amplitude as in The 
linear correlation in Fig. [21 suggests that an essentially R- 
independent result for the height of the FSI peak will be 
found if the LEC d is adjusted so as to realize a particular 
value of the Gamow- Teller matrix element in pp fusion. 
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FIG. 2: Gamow- Teller matrix element ||7|l plotted against 
TT~d nny FSI peak height given by JSJ for various val- 
ues of R. The nine points correspond to the values of R (in 
fm) indicated. The straight line is a linear fit to the points. 

This is implemented in xPT by adding to the matrix 
element Q the contribution from d: 
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Here the short-distance (r < R) pieces of the deuteron 
and pp radial wave functions are given by ^351 sin(K35ir) 
and AfgQ sin(K^ggr). We then extract a value for d by 
demanding that the total A^gt agrees with the result 
Ainj = 4.898 fm obtained from the Argonne Wig poten- 
tial [29j in Rcf. JJig 2M- (If choose the A^gt found 
using other NN potentials in Ref. the change in a„„ 
is < 0.01 fm.) This procedure gives the values for d{R) 
displayed in Table These values are "natural" |30| . 
Now d contributes to the FSI matrix element through 
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LEG 
FSI 



C— d A3S1K351 yli5oKi5o, (10) 
As 



where ^150 sin(Ki5o7') is the nn (radial) wave function 
for r < R. The TT~d — > nwy spectrum (both without 
and with the d contribution) for varying radii R is then 
plotted in Fig.l5l|3i|. 

Including d drastically reduces the dependence on the 
short- rang e physics that dominated the a„„ uncertainty 
in [13, |22| ■ If the entire spectrum is fitted our i?-related 
a„n uncertainty is reduced from ±1.0 fm to ±0.14 fm. 
Other, smaller, sources of theoretical error in the a„„ 
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FIG. 3: Calculated 7r~d — > nn'y neutron time-of-flight distri- 
bution, for two R values, without and with the LEC d contri- 
bution. The supplement of the 7n lab. angle is 9^ = 0.075 rad 
and Ann = —18.5 fm. The labels QF (quasi-free) and FSI in- 
dicate where the corresponding kinematics are dominant. 

TABLE I: Values of d for various matching radii R. Here we 
have used = 93.0 MeV and M = 939 MeV. 



R (fm) 


1.4 


1.8 


2.2 


2.6 


3.0 


d 


-1.27 


-0.559 


0.481 


2.07 


4.29 



extraction are: the role of P-waves in the nn final state, 
potential changes in wave functions when chiral two-pion- 
exchange is incorporated, and pion-range 0{Q'^) pieces of 
the two-body operator used to compute the spectrum of 
Fig. O 0, 123I . Now that it is clear how to remove the 
i?-related uncertainty wc will include these effects in our 
spectrum computation |25| . It should then be possible to 
extract a„„ from 7r~d nwy with a total theoretical un- 
certainty of < 0.3 fm if the entire spectrum is fitted and 
^ 0.05 fm if only the FSI peak is fitted. In the second 
case the theoretical error would be more than three times 
smaller than in any previous extraction [Illlllii. Ex- 
perimental uncertainties would be the dominant contri- 
bution to the a„„ error bar. Note that the values found 
for a„„ by fitting in both regions should be consistent 
within errors — a useful test of the theory. 

This demonstrates how xPT links apparently disparate 
reactions in few-nucleon systems. The correlations that 
result from this are important not just for their impact 
on the accuracy with which the nn scattering length 
is known. They also mean that measurements of, e.g., 
might help determine the strength of 
the chiral 3NF, and, if done with sufficient precision, 
could test the value of d extracted from tritium beta de- 
cay in 0. 

We appreciate discussions on these subjects with Mal- 
colm Butler, Kuniharu Kubodera, and Fred Myhrer. 
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